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We  present  transport  measurements  of  unshunted  dc  superconducting  quantum  interference  de¬ 
vices  (SQUIDs)  consisting  of  30  nm  wide  aluminum  nanobridges  of  varying  length  L  contacted  with 
two  and  three  dimensional  banks.  3D  nanobridge  SQUIDs  with  L  <  150  nm  (approximately  3-4 
times  the  superconducting  coherence  length)  exhibit  «  70%  critical  current  modulation  with  applied 
magnetic  field,  approaching  the  theoretical  limit  for  an  ideal  short  metallic  weak  link.  In  contrast, 

2D  nanobridge  SQUIDs  exhibit  significantly  lower  critical  current  modulation.  This  enhanced  non¬ 
linearity  makes  3D  nanobridge  Josephson  junctions  well  suited  to  optimize  sensitivity  in  weak  link 
SQUID  magnetometers  as  well  as  realize  ultra  low-noise  amplifiers  and  qubits. 


Josephson  junction  circuits  based  on  weak  links  are 
well  suited  for  applications  requiring  high  transparency, 
small  area  junctions.  For  example,  superconduct¬ 
ing  quantum  interference  devices  (SQUIDs)  constructed 
from  nanoscale  bridges  are  widely  used  in  magnetiza¬ 
tion  measurements  of  magnetic  crystals  and  ferromag¬ 
netic  clusters^— .  Furthermore,  weak  link  junctions  are 
attractive  candidates  for  coherent  quantum  circuits  since, 
by  construction,  they  are  free  of  lossy  tunnel  oxide — a 
potential  source  of  decoherence^.  However,  these  perfor¬ 
mance  gains  come  at  the  expense  of  reduced  nonlinearity 
in  the  current-phase  relation  (CPR)  when  compared  to 
conventional  tunnel  junctions.  The  CPR  in  weak  link 
devices  depends  strongly  on  the  junction  geometry^.  In 
particular,  thin  film  bridges  with  lateral  dimensions  com¬ 
parable  to  the  superconducting  coherence  length  £  are 
predicted  to  have  maximal  nonlinearity  when  contacted 
with  banks  of  greater  thickness  than  the  bridge  (3D),  as 
opposed  to  banks  of  equal  thickness  (2D).  The  3D  banks 
enhance  nonlinearity  by  acting  as  good  phase  reservoirs 
and  confining  the  variation  of  the  superconducting  order 
parameter  to  the  weak  link£. 

In  this  letter,  we  present  measurements  of  the  mag¬ 
netic  flux  response  of  unshunted  dc  SQUIDs  comprised 
of  aluminum  nanobridges  connected  to  2D  and  3D  banks. 
Our  measurements  show  enhanced  critical  current  mod¬ 
ulation  in  3D  nanoSQUIDs  and  confirm  that  thin  film 
aluminum  nanobridges  up  to  150  nm  in  length  («  3  —  4£) 
with  thick  contacts  approach  the  nonlinear  CPR  of  an 
ideal,  short,  metallic  weak  link.  Longer  devices  and  those 
with  2D  banks  exhibit  reduced  modulation,  in  quantita¬ 
tive  agreement  with  our  numerically  computed  CPRs. 
Thus,  3D  nanobridge  junctions  are  robust  nonlinear  ele¬ 
ments  for  use  in  nanoSQUID  magnetometers,  amplifiers, 
and  coherent  circuits  such  as  quantum  bits. 

Our  devices  consist  of  a  small  superconducting  loop 
interrupted  by  two  nanobridges.  Each  bridge  is  an  8  nm 
thick,  30  nm  wide  aluminum  film  of  variable  length  from 
75-400  nm.  These  dimensions  can  be  reliably  achieved  us¬ 
ing  electron  beam  lithography  and  a  liftoff  process.  For 
2D  devices,  a  single  electron-beam  evaporation  is  per¬ 
formed  at  normal  incidence  through  a  bilayer  shadow 
mask  to  simultaneously  define  a  narrow  bridge  and  wide 
banks  of  the  same  thickness.  To  create  3D  banks,  we  tilt 


FIG.  1:  (Color  online)  Atomic  force  micrographs  of  8  nm 
thick,  30  nm  wide  aluminum  nanobridges  contacted  with  (a) 
8  nm  thick  2D  banks  and  (b)  80  nm  thick  3D  banks.  Two 
nanobridge  junctions  are  incorporated  into  an  unshunted  dc 
SQUID,  as  shown  in  (b).  The  SQUID  loop  is  1.0  x  1.5  pm, 
and  is  calculated  to  have  5  pH  of  inductance. 


the  sample  in  situ  and  perform  a  second  evaporation  step 
to  deposit  80  nm  thick  banks  without  depositing  further 
material  in  the  narrow  bridge  region.  Atomic  force  mi¬ 
crographs  of  2D  and  3D  nanobridges  are  shown  in  Fig.  1. 
The  SQUID  loop  is  1.0  x  1.5  pm,  and  is  calculated  to  have 
Lioop  =  5  pH  of  geometric  inductance.  Each  chip  con¬ 
tained  six  simultaneously  fabricated  nanoSQUIDs  with 
bridges  of  varying  lengths.  The  devices  were  mounted 
in  a  copper  sample  box  equipped  with  a  superconducting 
flux-modulation  coil  and  measured  in  a  3He  cryostat  with 
a  base  temperature  of  265  mK.  All  wiring  ran  though 
discrete-element  and  distributed  lossy  filters  to  suppress 
spurious  electromagnetic  radiation.  Four-wire  transport 
measurements  on  all  six  SQUIDs  were  performed  using 
custom  electronics. 

In  the  inset  of  Fig.  2,  we  present  a  typical  current- 
voltage  (I-V)  characteristic  of  a  3D  nanoSQUID  with  75 
nm  long  bridges.  The  TV  characteristic  is  hysteretic, 
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FIG.  2:  (Color  online)  IcRn  product  as  a  function  of  bridge 
length  for  both  2D  and  3D  nanoSQUIDs  from  several  chips. 
The  dashed  line  in  the  main  figure  is  the  zero  temperature  pre¬ 
diction  for  a  short  metallic  weak  link  connected  to  ideal  phase 
reservoirs.  Inset:  a  typical  I-V  curve  for  a  3D  nanoSQUID, 
consisting  of  two  75  nm  bridges. 


and  the  retrapping  behavior  can  be  attributed  to  ther¬ 
mal  relaxation^.  We  observe  that  the  normal  state  re¬ 
sistance  Rn  increases  linearly  with  bridge  length,  with 
an  offset  of  «  1.5  for  3D  bridges  and  «  20Q  for  2D 
bridges;  these  numbers  are  consistent  with  the  normal 
state  resistances  of  the  banks.  We  subtract  this  off¬ 
set  and  plot  the  resulting  IcRn  product  in  Fig.  2  for 
both  2D  and  3D  nanoSQUIDs.  For  short  3D  devices 
(L  <  150  nm),  we  see  that  the  IcRn  product  saturates 
at  approximately  380/W.  For  comparable  2D  devices, 
the  data  are  clustered  around  250/iV.  For  longer  de¬ 
vices,  both  2D  and  3D,  we  see  a  monotonic  increase 
in  IcRn  with  bridge  length,  as  expected  for  a  long  su¬ 
perconducting  wire^.  To  interpret  the  data  for  short 
3D  bridges,  we  can  make  use  of  the  analytical  result 
from  Kulik  and  Omelyanchuk^  (KO-1),  who  calculate 
the  current-phase  relation  of  a  weak  link  with  L  -C  £ 
to  be  1(5)  =  (TtA/eRn)  cos(i5/2)  tanh_1[sin(<I/2)]  where 
5  is  the  phase  difference  across  the  junction  and  A  is 
the  superconducting  gap.  This  formula  assumes  a  diffu¬ 
sive  weak  link  at  zero  temperature  with  a  phase  vari¬ 
ation  constrained  to  the  bridge  region.  The  dashed 
line  in  Fig.  2  indicates  the  prediction  of  this  relation, 
IcRn  =  1.327rA/2e,  which  agrees  well  with  our  short 
3D  bridges.  For  our  thin  films,  we  estimate  the  coherence 
length  in  the  dirty  limit  to  be  £  =  sjKDj^nk^Tc  ~  40 
nm,  where  D  is  the  diffusion  constant  calculated  from 
the  measured  film  resistivity.  Thus,  our  shortest  devices 
are  a  few  times  the  coherence  length  and  approximate 
well  a  short  metallic  weak  link  connected  to  ideal  phase 
reservoirs. 

To  characterize  the  flux  dependence  of  the 


nanoSQUIDs  we  measured  I-V  curves  as  a  function 
of  applied  perpendicular  magnetic  field.  From  these 
data,  we  extracted  the  critical  current  Ic  as  a  function 
of  applied  flux,  and  have  plotted  this  quantity  for 
nanoSQUIDS  with  75,  150,  and  250  nm-long  bridges  in 
Fig.  3.  We  observe  that  the  critical  current  modulation 
is  deeper  and  more  nonlinear  for  the  3D  devices.  The 
modulation  for  the  2D  devices  is  quite  triangular  in 
shape,  indicative  of  a  nearly  linear  current-phase  re¬ 
lation.  The  data  show  a  modulation  pattern  which  is 
nearly  symmetric  about  the  maximum  Ic-  Comparing 
the  measured  modulation  with  numerical  simulations 
indicates  that  the  junction  critical  current  asymmetry  is 
typically  less  than  10%.  Larger  asymmetry  would  reduce 
the  SQUID  modulation  depth  and  thus  complicate 
systematic  comparison.  The  modulation  depth  in  our 
devices  is  also  not  limited  by  the  loop  inductance  as 
/3l  =  2LioopIc /‘ho  ~  0.1,  where  $o  is  the  flux  quantum. 
In  Fig.  4,  we  plot  the  fractional  modulation  depth 
( Ic,max  -  Ic,min)/Ic,max  as  a  function  of  device  length. 
For  short  3D  bridges  ( L  <  150  nm),  we  observe  a 
modulation  depth  approaching  70%.  In  longer  bridges, 
the  modulation  depth  decreases.  For  very  short  2D 
bridges  ( L  <  100  nm),  we  observe  a  modulation  depth 
below  50%.  In  intermediate  length  2D  bridges,  the 
modulation  depth  increases,  and  then  decreases  again 
for  bridges  longer  than  150  nm. 

To  interpret  these  results,  let  us  consider  a  SQUID 
formed  by  two  symmetric  elements  with  an  arbitrary  cur¬ 
rent  phase  relation  f(5).  In  the  limit  of  small  loop  in- 
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FIG.  3:  (Color  online)  Critical  current  as  a  function  of  applied 
magnetic  flux  for  75,  150,  and  250  nm  long  nanoSQUIDs.  The 
upper  panel  are  experimental  results  for  nanobridges  with  3D 
banks  and  the  lower  panel  for  2D  banks. 
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FIG.  4:  (Color  online)  Critical  current  modulation  depth  as 
a  function  of  device  length  for  both  3D  and  2D  nanobridges. 
The  black  dashed  line  is  the  theoretical  prediction  for  an  ideal 
weak  link  SQUID  at  zero  temperature  with  zero  loop  induc¬ 
tance.  The  blue  and  red  dashed  lines  are  numerical  calcula¬ 
tions  for  our  3D  and  2D  devices,  respectively,  including  the 
effect  of  finite  loop  inductance. 
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ductance,  the  phase  difference  between  the  two  SQUID 
junctions  52  —  equals  the  applied  flux  27r<E>m/<I>o.  The 
current  I  through  the  SQUID  is  the  sum  of  the  currents  in 
each  arm  I  =  /(<5i)  +  f(S2)  =  f(Si)  +  /(<5i  +  27r$m/<f>0)- 
For  a  tunnel  junction  based  SQUID,  f(S)  =  /osin(d); 
maximizing  the  current  equation  above  gives  the  SQUID 
critical  current  Ic  =  2/0|  cos(7r$m/<I>o)|  where  Iq  is 
the  critical  current  of  each  junction.  The  critical  cur¬ 
rent  modulates  from  0  to  2/o  and  hence  the  modula¬ 
tion  depth  is  100%.  For  a  linear  current  phase  rela¬ 
tion,  expected  for  a  long  wire,  the  maximum  modula¬ 
tion  depth  is  25%.  In  the  case  of  the  KO-1  relation 
cited  earlier,  the  maximum  modulation  depth  is  81%. 
Thus,  our  short  3D  bridges  approach  the  KO-1  result, 
while  the  longest  devices  tend  towards  25%  modulation 


as  the  CPR  becomes  increasingly  linear.  Also  plotted 
in  Fig.  4  are  the  modulation  depths  predicted  from 
our  numerically  computed  CPRs^.  These  predictions 
take  into  account  loop  inductance,  using  the  condition 
27 m  =  61  -  82  +  27T$m/4>0  +  TrLloop[f(51)  -  f(S2)]/$ 0. 
The  quantitative  agreement  with  short  2D  and  3D  de¬ 
vices  is  excellent,  except  for  the  shortest  2D  devices, 
which  modulate  less  than  predicted.  We  attribute  this 
to  the  fact  the  phase  gradient  in  short  2D  bridges  leaks 
into  the  banks  of  the  device.  Thus,  for  two  weak  links 
in  close  proximity  it  is  difficult  to  impose  a  phase  differ¬ 
ence  between  the  two  arms  of  the  SQUID&,  and  modula¬ 
tion  is  suppressed.  Our  calculations  are  for  isolated  junc¬ 
tions  and  a  numerical  calculation  of  the  full  two  junction 
SQUID  is  required  to  reproduce  this  effect^. 

In  summary,  using  low-inductance  symmetric  dc 
SQUIDs,  we  have  demonstrated  that  nanobridge  Joseph- 
son  junctions  contacted  with  3D  banks  exhibit  greater 
nonlinearity  than  those  with  2D  banks.  Transport  in  3D 
bridges  up  to  several  coherence  lengths  long  is  well  de¬ 
scribed  by  the  KO-1  formula  for  an  ideal,  short,  metallic 
weak  link.  Our  experimental  results  give  good  quanti¬ 
tative  agreement  with  a  numerical  solution  of  the  Us- 
adel  equations  for  the  precise  geometry  under  study.  As 
such,  3D  nanobridge  devices  should  have  sufficient  non¬ 
linearity  for  use  in  high  gain  amplifier  and  qubit  circuits 
traditionally  employing  tunnel  junctions.  These  devices 
could  also  be  used  for  dispersive,  high-speed  magnetom- 
etry  by  embedding  them  in  a  microwave  resonator—.  For 
magnetometry,  3D  nanoSQUIDs  should  have  higher  flux 
sensitivity  than  2D  devices,  but  perhaps  with  a  reduced 
tolerance  to  in-plane  magnetic  fields.  Further  study  is 
needed  to  determine  the  variation  of  nonlinearity  and 
field  tolerance  with  bank  thickness. 
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